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Abstract 

Dynamical Mean-Field Theory (DMFT) has established itself as a reliable and well-controlled 
approximation to study correlation effects in bulk solids and also two-dimensional systems. In 
combination with standard density-functional theory (DFT) it has been successfully applied to 
study materials in which localized electronic states play an important role. There are several evi- 
dences that for extended systems this DMFT+DFT approach is more accurate than the traditional 
DFT+U approximation, particularly because of its ability to take into account dynamical effects, 
such as the time-resolved double occupancy of the electronic orbitals. It was recently shown that 
this approach can also be successfully applied to study correlation effects in nanostructures. Here, 
we present a brief review of the recently proposed generalizations of the DFT+DMFT method. In 
particular, we discuss in details our recently proposed DFT+DMFT approach to study the mag- 
netic properties of nanosystemspQ and present its application to small (up to five atoms) Fe and 
FePt clusters. We demonstrate that being a mean-field approach, DMFT produces meaningful re- 
sults even for such small systems. We compare our results with those obtained using DFT+U and 
find that, as in the case of bulk systems, the latter approach tends to overestimate correlation ef- 
fects in nanostructures. Finally, we discuss possible ways to farther improve the nano-DFT+DMFT 
approximation and to extend its application to molecules and nanoparticles on substrates and to 
nonequilibrium phenomena. 

PACS numbers: 75.10.-b, 75.75.-c, 75.40.Mg, 75.50.Bb, 75.50.Cc 
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I. INTRODUCTION 



The theoretical description of solid state systems which include localized electronic states 
(usually in the case of transition metal (TM) atoms with unfilled shells of d-electrons) re- 
mains one of the most difficult problems of condensed matter physics. These states play an 
important role in many unusual phenomena in different materials, like cuprate supercon- 
ductors, heavy fermions, manganites and others. In the case of small TM systems, which in 
modern condensed matter language correspond to nanosystems, clusters and molecules, the 
role of electron-electron correlations are of special interest from both practical and fundamen- 
tal points of view. Indeed, these structures have a great potential to be used in many modern 
technological applications, like ultra-high density magnetic recording systems, sensors and 
photovoltaic devices. From the fundamental point of view, understanding of correlation 
effects in small TM systems is important because it will help understand growth dynamics 
and the role of the quantum effects in the corresponding larger structures. Also, the lack of 
experimental data for many such small systems places a great burden on theory for predict- 
ing their properties. Unfortunately, most studies of TM clusters either ignore strong electron 
correlation effects, or take them into account in an oversimplified way, thereby neglecting 
the very effects that maybe responsible for their unusual properties such as magnetism and 
superconductivity. In the case of systems at the nanoscale correlation effects are expected to 
be even more important due to reduced dimensionality, which makes it much more difficult 
for electrons "to avoid" each other as compared to the case of bulk material. 

Traditionally, in the case of nanostructures and molecules correlation effects are included 
using the DFT+U approachj2j |3]. However, in this (mean-field) approximation important 
dynamical effects such as time-resolved local interactions are neglected. These effects appear 
to be rather important when correlations are not extremely large (small on-site Coulomb 
repulsion energy U), as in the case of bulk TMs, for which the DFT+U method may lead 
to wrong results, such as predicting a non-existing magnetic phase in plutonium. In the 
case of large U, the DFT+U approximation may also lead to wrong results, for example 
spin-ordering temperatures (see, e.g., review Ref. [I] and references therein). 

On the other hand, in the case of extended systems several many-body theory approaches 
have proven to incorporate electron correlation effects successfully. Familiar examples in- 
clude the Bethe ansatz approach in the ID case (see, e.g., Ref. [3]) and the DMFT method 
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for higher dimensions. In the last approach, one takes into account temporal fluctuations 
by considering the time- (frequency-) dependence of the electron self-energy. The spatial 
dependence of the self-energy is neglected, which strictly speaking is exact only in the limit 
of infinite dimensions (coordination number Z)[6J. Nevertheless, DMFT appears to be a 
good approximation for 2D and 3D systems (for an over-review, see Ref. [7j). The DMFT 
formalism has also been successfully extended to incorporate electron-phonon interaction [8] , 
nonhomogeneous systems (multilayered structures) [9HTTj. and nonequilibrium phenomena 

An important step in the analysis of the correlation effects in real materials was made 
when it was proposed that the DMFT approach can be incorporated into the DFT scheme 
[TSUll]. In this combined DFT+DMFT method, the system geometric structure (lattice, in- 
teratomic distances etc) and the properties of corresponding "non-correlated" (Kohn-Sham) 
electron subsystem are obtained within the DFT approximations (usually LDA or GGA), 
and correlation effects are taken into account via the short-range Coulomb interaction of the 
quasiparticles within a Hubbard-type tight-binding model. In the case of bulk materials, 
this approach has been successfully applied for studying the spectral, optical and magnetic 
properties of systems (see Refs.jl] and [15] for reviews). In particular, it was shown that 
the dynamical fluctuation effects which are naturally incorporated in DFT+DMFT approach 
may lead to a different behavior of the system from that arising from DFT+U, as mentioned 
above. 

At first sight, it is not obvious that the DMFT approach can also be applied to nanos- 
tructures or molecules, for which space nonhomogeneity is important in most cases. It has 
been suggested by Florens [T6] . however, that when the average atomic coordination number 
(Z) is large, DMFT is a good approximation even in the case of finite-sized systems. Tha 
application of DMFT to nanostructures, namely, surfaces, nanoparticles, molecules etc, may 
also be feasible since Z in these systems ranges from vary small to bulk values. The main 
contribution to the spatial fluctuations comes from the difference between the "surface" 
(low-coordinated) and the "bulk" atoms. Therefore, one can expect that the spatial fluctu- 
ations are important when the number of the surface atoms is comparable to or larger than 
that of the bulk atoms. This is definitely not the case for nanostructures (except chains), 
so the DMFT approach must be a reasonable approximation for these systems. Moreover, 
it is expected to be a good approximation even in the case of much smaller systems, down 
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to non-flat 5 — 10 atom clusters, since in this case the average coordination number is of 
order of 4, similar to that one of the 2D systems, where DMFT is often considered to be 
a good approximation. Furthermore, there is the generalization of DMFT for the nonho- 
mogeneous case, in which the electron self-energy remains local, but site-dependent j9j ITU] . 
This generalization and enhanced accuracy, which is computationally much more time- and 
memory-consuming, allows one to study systems with several hundred atoms within available 
computational resources (see Section II). 

The DMFT and combined DFT+DMFT approximations were recently applied to study 
several properties of systems containing up to one hundred atoms [U [T7H2T] . For example, 
spectral properties of the Mn4 and ring molecules were analyzed by using DFT+DMFT. 
[TT] 120] It was shown that the DMFT approximation not only gives resonable results even 
for such small systems, but also a better description as compared to the DFT+U [17] or 
unrestricted HF[20J approximations. The transport properties of 110 atom quantum point 
contact [IB] and small Ni clusters between Cu nanowires [TH] were analyzed by using DMFT 
and DFT+DMFT, respectively. For quantum chemistry, as an alternative to DFT+DMFT 
a combined Hartree-Fock (HF) and DMFT approximation was used to study the spectral 
properties of solid bulk hydrogen. [21 J The effects of the short range interaction were ex- 
amined within DMFT, while the remaining long-range interactions were treated with the 
HF approach along these directions. We have recently proposed a combined DFT+DMFT 
approach to study the magnetic properties of nanosystems, and showed its validity for ex- 
amining the magnetic properties of small iron clusters [T]. 

In this paper, we provide details of our recently proposed combined DFT-DMFT approach 
to study the physical properties of finite systems, which include nanoclusters and molecules, 
and in which electron-electron correlation play an important role. We elaborate on the 
application of the method for nanomagnetism and provide results for the magnetic properties 
of very small (2-5 atoms) Fe and FePt clusters, We analyze the role of correlation effects for 
different cluster geometries and chemical composition. Our comparison of the DFT+DMFT 
and the corresponding DFT+U results suggests that neglect of dynamical correlations tends 
to lead to an over-estimation of the role of electron correlations. Finally, we discuss possible 
generalization of the approach to nanostructures and molecules on substrates and to non- 
equilibrium systems. 
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II. THE DMFT APPROACH: FROM BULK TO NANOSIZE 



In this Section, for completeness we first review the main equations for the DMFT formal- 
ism for extended (homogeneous) systems and its (nonhomogeneous) generalization for finite 
systems. In most cases, the properties of systems with strong electron-electron correlations 
can be studied by solving the problem with the tight-binding Hamiltonian of the Hubbard 
type: 

H ^ ^ ^iiy'mCj^jCjcrm ^ U a{r iTli (T iTlj a i rn . (1) 

i,j,cr,l,m i,j,cr,a',m 

In the last expression c\ al and c Jcrm are the electron creation and annihilation operators 
of the fermion at site i with quantum numbers o (spin) and and 1 (other quantum numbers: 
orbital momentum, band number etc). In these notations, the particle number operator, 
which corresponds to state a, / is n ia i = c\ al c i(T i. The competition between the kinetic energy 
and the (short-range Coulomb repulsion) potential energy is defined by the hopping t^ym an d 
the Coulomb repulsion matrices U l ^ aa i. Usually, the hopping matrix elements are obtained 
from fitting the experimental band structure, and only the nearest-neighbor or next-nearest- 
neighbor hopping matrix elements are taken into account. The Coulomb repulsion energy 
is often chosen in the local approximation, (U^ aa , = 5ij5i m U l ), since these terms give the 
largest contribution into the Hamiltonian in Eq.Q, and the values for U l (/ is usually the 
orbital number: s, p, d ...) are chosen as fitting parameters, though one can in principle 
obtain the values for U l from DFT calculations (see, e.g., Ref. [15]). 

The problem described by the Hamiltonian in Eq. Q can be solved by using the time- 
ordered Green's function (GF) 

(Wm(M') = -i(Tc itTl (t)cl m (t')). (2) 

This function is connected to the single-particle self-energy Tii a i-j a ' m (t, t') through the Dyson 
equation in the following way: 

G^Jyj) = GfLusJu) + [G^\u)E(u;)G(u;)]. al;jalm (3) 

(in the frequency representation). In the last equation, Gf^.-^^u) is the non-interacting 
GF (which corresponds to the {7 = case). Since in the high- dimensional (coordination) 
limit, the problem can be reduced to the single-site problem with the local GF and the local 



self-energy, taken for example at sites i=j=0, one can find from the Dyson equation ^ the 
following expression for the local GF: 

G aWm (u) = [ ( 77Y7 ^TT ) • ( 4 ) 

J (27r) rf \u - e(k) + n - ) al . a , m 

In the last equation, e(k) is the "free" quasi-particle spectrum obtained from the Fourier 
transform of the hopping matrix, \x is the chemical potential, and is the quasi-particle 
self-energy, which is also local and site-independent in the case of the DMFT approximation. 
The frequency-dependence of the self-energy allows one to take into account dynamical 
effects, such as the time-dependent local occupation numbers and interaction. In the case 
of standard mean-field (in particular HF, which corresponds to DFT+U) approximation, 
in which S is assumed to be static (frequency-independent), these effects are neglected. 
Though finding the frequency-dependence of the self energy is often not an easy task, the 
space- (momentum-) independence of the self-energy is a crucial simplification, which is exact 
in the limit of infinite dimensions or coordination number [6], and allows one to solve the 
problem of strongly correlated electrons with rather high accuracy in many cases. A variety 
of physical quantities can be obtained from the GF, for example the orbital spin density 

" ri = ~/^/ J^ lmG ^ u) - (5) 
It follows from Eq. Q that in order to find the GF one needs to find the self-energy 
E(tu). Additional equation(s) which connect the local GF and self-energy can be found in 
the following way. The local GF can be formally expressed in terms of the path integral 
over the Grassmann variables ijj and ijj* : 



<Wm(w) = J D[^]D[r]M:> m e- A ^ r \ (6) 
where the action A has the following form: 

(here and in some places below, for simplicity we neglect the orbital degrees of freedom and 
assume the Coulomb repulsion to be local). The intergation over all states which correspond 
to the sites different from I = in Eq. ^ can be formally performed, which gives: 

G^Ju) = [ Dty]D[r}Air*> m e- Aeff[ ^ r ' S ~ 1] , (8) 



where 

A eff W,p,g- 1 ] = - I dr [ dr'^MS-'ir-r'^oair^ + U I drn^{r)n^{r) (9) 
JO Jo a Jo 

is the impurity effective action (r is the imaginary time and /3 = 1/T is the inverse temper- 
ature). This action depends on the effective dynamical mean-field Gai-ymiu), which takes 
into account all effects of the rest of the system with sites ^ on the impurity site. The 
problem now is equivalent to the problem of the single site coupled to the bath described 
by the field Go-i;o-> m (u)). In principle, in the case of given dynamical mean- field one can find 
the local GF from Eq. ([8| by using different approaches (see below). 

The next step is to find the dynamical mean-field Goi;a'm{w)- The equation which connects 
this field with the local GF and self-energy can be found by mapping the single impurity 
problem to that of the Anderson impurity with the Hamiltonian: 

h a = 5^(eo - n)cl a c 0(T + U n ot n ol + ^ V k (c ] K(T c 0(7 + cJ^Ck^J + ^ 44,^1^- (10) 

cr k,tr k,(T 

where eo and e k are the energy levels of the impurity and of the bath electrons, and Vk 
is the hybridization between the impurity and bath states (see, e.g., Ref. [7]). Since this 
Hamiltonian is quadratic in the bath-field operators, one can easily find the impurity Green's 
function in the non-correlated case (U=0): 



where 



a m = £tS> ( 12 ) 



, u - e k 

k,tr K 



is the hybridizarton function. Since, we assume that "the Hubbard" (in Eq. Q) and "the 
Anderson" (Eq. ( 10 )) baths are equivalent, the mean-field functions in Eq. ^ and in Eq. ( 12 ) 
are equal to each other (since in both cases they are equal to the local impurity GFs at 



U=0). On the other hand, the function (12) is connected with the impurity self-energy in 
the standard way by means of the Dyson equation 

G- l {u) = g- l {u)-Y,{uj). (13) 

Again, since the two impurity problems are equivalent, the self-energies in Eqs. Q and ^ 
have to be the same. 
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We thus obtain a closed system of DMFT equations, i.e. Eqs. Q, ([8| and (13), for 
the local GF, dynamical mean-field and self-energy. It can be solved by iterations for every 
value of frequency in the following way: 

1. choose the initial self energy E(cu); 

2. calculate the local GF from Eq. Q; 



3. calculate the dynamical mean-field G(oo) from Eq. (13); 

4. use the mean-field function to solve the impurity problem ^ to find the local GF; 

5. find new self-energy from Eq. (13) by using the values of G(u) and G(uj); 



6. continue iterations until S(cu) converges with the desired accuracy. 

The most difficult part of the solution is the path integral equation ([8]). Usually, exact 
numerical methods like Quantum Monte Carlo (QMC) solvers (e.g., Hirsh-Fye or continuous- 
time QMC approaches), or analytical ones employing expansion either in powers of t/U or 
U/t, including the iterative perturbation theory expansion (for details, see Refs. [H [71 ITo]). 
are used. 

In the case of nanostructures, the problem is essentially nonhomogeneous (except some 
specific cases, like symmetric clusters that consist of single type of atoms: dimers, equilateral 
triangles etc., in which all sites are equivalent), making the process complex and tedious. To 
simplify the solution, it was proposed [16] that the properties at a given site (a cavity) can 
be obtained by knowing the solution on the rest of the sites. To find the last solution, one 
can continue by considering the next cavity approximation, by removing the nearest to the 
cavity shell and so on. Thus, when one finds the solution for the external shell one can move 
backwards, to solve one-by-one the internal shell problems coupled to the corresponding 
external shell bath. At the end, one can find the dynamics on the central atom (see Fig. 1, 
where the illustration of the reduction of the problem to the cavity problem for a three-shell 
nanoparticle is presented). 

In a more straightforward manner, in the case of not too large system (up to a few 
hundred of atoms) one can solve the problem "exactly" without dividing the system into 
shells by generalizing the system of the DMFT equations in real space representation (the 
generalization of the DMFT theory on the nonhomogeneous Hubbard model was made by 
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FIG. 1: The illustration of the reduction of the DMFT problem to the cavity problem in the case 
of three-shell nanoparticle. 

Nolting and Potthoff [9l [in])- Namely, the local site-dependent GF for the N-atom system 
can be obtained from equation (|3]), which can be written as: 

+ t 12 ti 3 ... t 1N \ 

tlx W + /i-S 2 (w) t 23 ••• t 2N 

Cr , (.c)= t 31 t 32 oj + h-E 3 (uj) ... t 3N ■ (14) 



\ t N1 t N2 t N3 ... W + \L - Ejv(w) / 

We assume that the self-energy remains local in space, but it is site-dependent: T^jioj) 



5ijY,i(u). The other two equations (Eqs. ([8]) and (13)) remain the same, as in the extended 
case, though one needs to solve N impurity problems (or less, if some sites are equivalent 
due to the symmetry of the system): 



(15) 



G^) = J D[^\D[r\^*^v(- J^dr jf^^CW^fr-^M 



+U J drn it (r)fi4(r) ) . (16) 



Equation (16) can be obtained in the same way as in the bulk case, when instead of one site 



one considers every site of the system: 1 < i < N. Equation (15) is less obvious. To derive 
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it, one can map the single impurity problem to the Anderson impurity problem for every 
site of the system. In this case, one obtains N bath functions Gu(uj). 

Therefore, in the case of nanosystems one needs to solve the system of DMFT equations 



(14), (15), and (16). Similar to the extended case, one can solve the problem by iterative 



procedure for every value of frequency: 

1. choose the initial self-energies £j(u;) for each site; 



2. find the local GF, Ga(u), from Eq. (14) for each site by inverting the matrix and 



extracting the diagonal elements (zth diagonal element is equal to Gu(oj)); 



3. calculate the mean-fields Qu(u), 1 < % < N from Eq. (15) 



4. use the mean-field functions to solve the impurity problems (16) to find the local GFs 
for all sites; 



5. find new self-energy from Eq. (15) by using the values of Gu(oj) and Gu(oj); 



6. continue iterations until £^(u;), 1 < i < N, is converged with the desired accuracy. 

In both extended and nano-cases, the problem can be solved separately for every fre- 
quency, which makes it easy to parallelize. However, the nonhomogeneity in the nano-case 
makes the calculations far more demanding as one needs to work with the inversion of rather 



large matrices (14). Assuming that modern computers are capable of inverting efficiently 
matrices of size up to 10.000 x 10.000, systems with 500 atoms or less can be studied by using 
"the exact" set of equations. This number is obtained by dividing 10.000 by the numbers 



of orbitals and spin projections which also must be included in the matrix (14). Usually, 
one takes around 6 or 9 orbitals into account, which may include one s-orbital, three p- 
orbitals, and five d-orbitals etc. This allows one to study nanoparticles of size up to ~ 5nm. 
For larger systems, one could divide them into blocks, solve for each block separately, and 
then sew the solution by solving the problem of the "border" layers between each block, by 
"freezing" the internal (core) block sites, i.e. by keeping the local self-energy for these atoms 
fixed. The accuracy of the DMFT solution can be estimated by comparing the results with 
some exact results known for the Hubbard model. In particular, one can compare the values 
of the lowest Green's function and self-energy spectral moments with the exact analytical 
expressions. [22TT24] 
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III. DFT+DMFT APPROACH: BULK CASE 



In the case of DFT+DMFT calculations the first step consists of a DFT study of the 
system, i.e one needs to optimize the system structure and obtain characteristics of the 
"non-correlated" quasi-particles, such as the bandstructure (energy e(k) in Eq. Q)) and 
the density of states (DOS). At the next step, one can extract the parameters for the tight- 
binding Hamiltonian: the inter-site hopping parameters and short-range (usually, local and 
nearest-neighbor) Coulomb repulsion energies. Next, the correlated problem can be solved 
allowing estimation of the role of correlations in properties of interest. To make the analysis 
more consistent one should perform the DFT calculations by using the DFT+U method, to 
ensure that at least the non-dynamical part of the correlation effects is taken into account 
when optimizing the system geometry and that the band structure is obtained with the 
same value of U that is used in the DMFT calculations. However, one must be careful 
not to take into account the correlation effects twice: in geometry optimization and in the 
DMFT calculations. The optimization process is not affected by this double counting, since 
it is performed only one time, but the band structure in principle must be corrected by 
extracting the contribution which comes from the LDA or GGA exchange-correlation (XC) 
potentials. In many cases this contribution is not very large and may be neglected. 

The DFT+DMFT approach continues to be successfully applied to many correlated elec- 
tron materials. We refer the reader to excellent reviews of this topic [U [T3] , and to more 
recent numerous literature, most notably that of iron pnictide superconductors. [251 

IV. DFT+DMFT APPROACH: NANOSYSTEMS AND MOLECULES 

In the case of small systems, the DFT+DMFT formalism has the same strategy as in the 
bulk case: first, to obtain the main properties of "non-correlated" electrons within DFT, or 
more precisely within DFT+U, and then to use the resulting structures and other quantities 
to construct and solve the tight-binding Hamiltonian. 

Contrary to the bulk, in the nonhomogeneous/nanocase it is more convenient to solve the 
problem in real space. Therefore, one needs to extract the effective hopping parameters from 
the DFT calculation. In general, the hopping parameters are equal to the matrix elements 
of the non-interacting Kohn-Sham Hamiltionian with respect to the localized atomic-orbital 



11 



wave functions: 

r v 2 

Ujafi = / drip* a (r + Ri)(-— + V atfmlic (r))^(r + R,), (17) 

where R« and Rj are the site vectors and a and /3 are the corresponding orbital numbers, 
and V atom i C (r) is the atomic potential. Similar, to the bulk case, the values of the Coulomb 
repulsion parameter can be obtained either from the DFT calculations or it can be used as 
a parameter. 

Below, before giving a relatively detailed description of our application of this approach to 
magnetism (Subsection IV. B), we briefly summarize other applications of the DFT+DMFT 
approach to study the spectral and transport properties of nanosystems and molecules. 



A. Spectral and transport properties 

Perhaps for the first time the electronic structure of small systems was studied within 
a DFT(LDA)+DMFT by Boukhvalov et al HZ], who used the cluster DMFT approach to 
study properties of Mn^ molecular magnet. It was shown that DMFT predicts correct 
electronic gap observed experimentally in optical conductivity measurements. The LSDA 
and LDA+U calculations result in a finite DOS at the Fermi energy and a much smaller gap 
(0.9eV), correspondingly (in the LDA+U case the same values of Coulomb repulsion and 
the intra-atomic Hund exchange, U=4eV and J=0.9eV, were used). As these calculations 
show, the dynamical correlation effects can play a crucial role in the electronic structure of 
small systems. 

The cluster DMFT was also used to analyze the ground state geometry and the on-site 
projected DOS of the hydrogen clusters (H4 cluster, if 6 chain and ring, H 50 chain) [20J. 
As the authors have shown, in the case of intermediate values of Coulomb repulsion the 
DMFT is a superior approach comparing to other methods, including the unrestricted HF. 
In particular, DMFT predicts reasonably accurately the correct position of the minimum of 
the energy of tetragonally-coordinated H$ cluster as a function of the inter-atomic distance, 
while the unrestricted HF approach fails in this case. 

The system energy and the DOS of the extended cubic system of hydrogen atoms was 
studied in Ref. [2T]. In particular, the dependence of the properties of the system on the 
number of bath orbitals was analyzed, which is an important question in the case of finite 
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systems. 

In Ref. [TH], the diagrammatic dynamical approach to solve the Hubbard model coupled to 
a non-interacting environment was used with applications to study the transport properties 
of different nanosystems, like a few-atom system, benzene molecule and a quantum point 
contact consisting of 110 atoms. In particular the spectral properties and conductivity of 
the systems was analyzed by taking into account the possibility of the Mott transition due 
to the electron correlations. This approach was proposed as an alternative to avoid using 
the Cayley-type of the model, discussed by Florens[IB], where the properties of the central 
site depend on the properties of the external ones without inverse dependence. Notably, 
as the authors have shown the n-vertex approximation gives meaningful results already at 
n=l, i.e. at the DMFT level. 

The transport properties of nano-contacts were studied in Ref. [19J, where the DMFT 
GFs were used to calculate the Landauer transmission function. Namely, the case of Ni 
nanocontacts between Cu nanowires was analyzed and it was shown that the dynamical 
correlation effects can significantly alter the behavior of the system. 

B. Magnetism 

To solve the nano-DMFT problem for specific Caley-tree-type system, one may use Flo- 
rens algorithm. On the other hand, for small systems (up to few hundred atoms) we can 
solve the problem exactly, as outlined above. For systems containing few atoms, one can 
also use DMFT codes for extended (periodic) systems, by assuming, for example, that every 
cluster occupies a site on a two-dimensional lattice with a very large lattice constant such 
that the cluster-cluster interaction may be neglected (a super-cell approximation). In this 
case, every atom in the cluster corresponds to the site orbital (which includes the atomic 
orbitals as an additional quantum number) . We have applied such an approach to study the 
magnetic properties of small Fe and FePt clusters. The DMFT part of the calculations was 
performed using the LISA code [I] with the Hirsch-Fye QMC solver [27] . These calculations 
were based on the cluster structures obtained by a DFT calculations (VASP4.6 code[28j, 
GGA-PW91 XC potential, for details, see Ref. [T]). To take correlation effects into account 
starting at the DFT level, the DFT+U approach was used, where we assumed that U is site- 
and orbital-independent (3d- and 4s-states were taken into account for Fe atoms, and 5d- 
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and 6s-states in the case of Pt atom). The examples of the cluster geometries used in the 
DMFT calculations are presented in Fig. 2. As reference value for the Coulomb repulsion we 
have used U=2.3eV, often employed in the case of bulk Fe. However, the properties of the 
system were analyzed by varying U from to approximately 5 eV. In particular, to check 
our DFT+U results with experimental data, we compared the bondlength for the Fe dimer 
for U=2.3 eV with the experimental value, the only available data for small Fe clusters. For 
this case, our calculations find bondlength of 1.99A, with DFT+U in good agreement with 
the experimental estimations 1.87 A [2H] and 2.02A [50] . 

From the DFT+U calculations we have obtained values for the hopping parameters for the 
s- and d- valence electrons by using the Slater-Koster matrix approximation [311 132] . in which 
the effective d-orbital radius was approximated by the Muffin-Tin-Orbital values 0.864A 
and 1.1 16 A for the Fe and Pt atoms, respectively (see, e.g., Ref. [32]). In the case of iron 
clusters, we rescaled the hopping parameters by a factor of 0.367 in order to reproduce the 
average experimental value for the spin- up and spin-down bandwidth, which are split, in this 
case. This rescaling may also be considered equivalent to changing U. Indeed, it is not the 
individual values, rather the ratio between the hopping parameters for the different orbitals 
for the atoms in question that needs to be kept constant. This also follows from the fact that 
all the quantities depend on U/r Z J 2 and U/r^, as it follows from the renormalization of the 
Hamiltonian. To make the calculations consistent, we have used the same renormalization 
also for the FePt clusters. The following periodic (cluster) DMFT calculations for the 
optimized structures were performed by allowing a weak inter-site (inter-cluster) hopping, 
for faster convergence. For this purpose, a weak staggered external magnetic field was also 
used. The calculations appear to become slow at large values of U and low temperatures. 
However, to obtain system properties in the zero electronic temperature limit, one does not 
need to necessarily consider the case T=0. When T is much lower than the hopping and 
Coulomb repulsion parameters - the other energy scales of the system - the properties of the 
system will be very close to the zero-temperature case. 

C. Application of nanoDMFT to small Fe and FePt clusters 

The magnetism of small Fe clusters has been studied by using several approaches, includ- 
ing DFT (see, e.g., Ref. [33]) and DFT+U (see, for example, Refs. [34TI36] ) approximations. 
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TABLE I: The magnetization per atom as function of the number of atoms in the Fe clusters at 
U=2.3eV. Along with the DFT+U and DFT+DMFT results we present a set of experimental data 



N DFT+U DMFT exp. 

2 4.00 2.03 3.25 ±0.5 

3 3.33 2.22 2.7 ±0.33 

4 4.00 1.84 2.7 ±0.8 

5 3.60 1.98 3.16 ±0.33 



Unfortunately, the experimental data on such small clusters is rather limited[29l [30J, 13?] . 
though it is sufficient for estimation of the validity of the DMFT approximation. In this 
Subsection, we first summarize the results of our analysis of the magnetic properties of 
small Fe clusters, part of which was presented in Ref. PQ, and then focus on the properties 
of bimetallic FePt tetramer. 

In the Fe case, we have considered clusters containing 2 to 5 atoms and several possible 
geometries, which correspond to the minimum energy configurations at different values of 
U. In particular, at U=2.3eV we obtained bondlength of 1.99A in the case of the dimer (the 
bulk iron bondlength is 2.49A). Our minimum energy structures also include an equilateral 
triangle, a single-side pyramid for four atoms and a bi-pyramid for five atoms (Figure 2). The 
calculations of the magnetization with the DFT±U approach show that the magnetization 
per atom is larger than the bulk value 2.2/xb- AS can be seen from the results plotted 
in Fig. 3, the magnetization grows as U increases and becomes U-independent after some 
critical value of the Coulomb repulsion ~ leV, which depends non-trivially on the number 
of atoms in the cluster. 

In Figure 4, comparison of the results for Fe^ and Fe^ obtained with both DFT±U and 
DFT±DMFT methods are presented, the behaviour is qualitatively the same for clusters of 
other sizes. As it is evident from this Figure, dynamical effects play a role even at rather 
small values of U, leading to a significant reduction of the cluster magnetization. 

In Table I we present results for cluster magnetization as a function of the number 
of atoms at characteristic value of U=2.3eV (for more details see Ref. [Tj). Once again, 
the DFT±U calculations tend to over-estimate cluster magnetic moments and dynamical 
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FIG. 2: Examples of the structures of the Fe (top row, U=2.3eV) and FePt (middle and bottom 
row, U=4eV) clusters obtained from the DFT+U calculations, that were studied in details in the 
paper. The red and the gray balls are the Fe and the Pt atoms, correspondingly. The bondlengths 
are given in Angsroms. In the case of bimetallic clusters we have presented the lowest energy 
geometries for both 3D (middle row) and planar (bottom raw) cases. 

effects in general lead to a reduction of the magnetization. While neither method gives 
perfect agreement with experiment, DFT+U consistently overestimates and DFT+DMFT 
somewhat underestimates the magnetization. The point here was not to obtain perfect 
agreement with experiment, rather to show that even for very small clusters, DFT+DMFT 
gives reasonable results. One would naturally expect that the DMFT approach will be much 
more accurate in the case of larger clusters, which would make it preferable over DFT+U. It 
must be mentioned that the values of magnetization for the selected small clusters used here 
as a test of our method are also reasonably well reproduced by other DFT approaches j38l-fld~] . 
and most calculated values lie within the experimental error bars, which are quite large. [37] 
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FIG. 3: DFT+U results for the magnetization as a fuction of U for the Fe2 — Fe§ clusters. 

Finally, for the small Fe clusters we would like to mention briefly our analysis of the the 
role of the dynamical effects in the case of Jahn- Teller distorted clusters. As our calculations 
in the case of Fe trimers show, the dependence of the magnetic moment on the cluster geom- 
etry is qualitatively the same in both DFT+U and DMFT cases, but similar to the previous 
results DMFT leads to a reduction of the cluster magnetic moment pp. In general such 
a significant reduction of the cluster magnetization may indicate that the orbital position 
and/or their occupancy may change dramatically (through the electron self-energy) when 
dynamical effects are taken into account. 

Theoretical results for the magnetic properties of small FePt clusters is even more limited 
than those for the case of Fe. (121 H3] Since there is hardly any experimental data available, 
it is very important to perform a systematic ab initio analysis of such small structures, in 
particular the evolution of their properties with size in order to understand the structure- 
function relationship. These clusters are of special interests for several reasons. For example, 
there is the need to understand how the addition of the higher, 5d, localized states of the 
nonmagnetic metal element affects the magentic properties of the 3d states of the magnetic 
atoms in the cluster. We have chosen Fe- and Pt-atom clusters as our test cases, since such 
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nanoalloys have been the subject of both technological and fundamental interest. 

In the case of the FePt clusters, we focused on the dimer and the tetramer. The first 
case allows one to analyze qualitatively the role of the interplay of the 3d and 5d orbitals 
in cluster magnetization. The second case was examined to understand the dependence of 
magnetism on the chemical composition of the cluster. Similar to the Fe case, the clusters 
were optimized by using the DFT+U approach, but contrary to the iron case it was found 
that the geometry of the clusters strongly depends on the value of U (Figure 2). In some 
cases two-dimensional (2D) structures have lower energy. In particular, the relative energies 
for the corresponding structures with the same values of Fe and Pt atoms at U=4eV are: 
FeiPh - planar: -0.13eV, 3D: OeV; Fe 2 Pt 2 - planar: -0.63eV, 3D: OeV; Fe 3 Ph - planar: 
OeV, 3D: -0.13eV. The average bond length is longer in the 2D geometry as compared to 
the 3D ones, which would make correlation effects more important for them. This follows 
from the notion that one may assume that in the 2D case less orbitals are involved in 
the hybridization, and thus correlation effects are enhanced because of the contribution of 
unbonded charge. 

Since the cluster geometry and the energy levels used in the DMFT analysis are obtained 
with the DFT+U calculations, we first analyzed how these quantities are affected by chang- 
ing U. The DFT+U results for the magnetization of Fe and FePt dimers as function of U 
presented in Figure 5, show that the magnetization barely depends on U for the latter (the 
value of the magnetization per atom is in agreement with other DFT calculations |43] ) . 
but it has a sharp increase at U ~ 3eV, in the case of the former. This result can be 
explained as follows. Since increase of U leads to an increase of the bondlength, and this 
value of U is of order 3eV - a typical value which corresponds to the Mott transition in 
bulk materials, and the change of the cluster magnetization is rather sharp at this point, 
one may speculate the presence of some kind of (magnetic) transition. The fact that the 
magnetization is U-independent in the case of FePt cluster can be explained by assuming 
that the d-orbitals of the FePt are much less hybridized in the last case as compared to the 
Fe dimer due to the energy mismatch even at small U, so increasing U which leads to an 
increase of the bondlength does not lead to a significant charge redistribution (localization) 
as compared to the case of the iron dimer. In fact, our DFT+U calculations give the fol- 
lowing results for the bondlength (at U=2.3eV): 1.99A (Fe 2 ) and 2.29A (FePt). In the case 
of increasing U the bondlength increases in both cases, for example in the case of FePt the 
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bondlength is equal to 2.34A at U=4.0eV (In the case U=0eV, the DFT calculations |43j 
give the bondlength 2.2 A). 

To get a better understanding of the role of the correlation effects in the spectral properties 
of the dimers we have calculated the positions of the s- d-state peaks as functions of U 
for both clusters (Fig.6). We find that with U increasing the distance between the spin- 
down peaks (two red lines) and spin-up d-orbital peaks (top blue solid line) significantly 
increases when U changes from 2eV to 3eV in the case of iron dimer. This implies that 
the the probability for one electron to migrate back and fourth from the spin up to the 
spin down orbital, and hence to decrease the cluster magnetization, decreases at large Us. 
This transition may correspond to the magnetization jump in Fig.5. In the case of FePt 
cluster, there is no such sharp change between the energy difference for these states (the 
corresponding dashed lines). Therefore, the magnetization does not change significantly with 
U. There is another interesting effect of changing U in the case of FePt cluster. Namely, the 
positions of the spin up d-state and spin-down s-state change when U is between 2eV and 
3eV. This means that at large Us one can neglect to some extend the role of s-orbital in 
correlation effects for FePt cluster, since this orbital may be assumed to be doubly occupied 
(this statement is true only when the distance between the s-spin-down and d-spin-up peaks 
is large enough, though). 

To obtain a deeper understanding of the role of different d-orbitals in the dynamical 
processes, we have analyzed the positions of projected d-states at different values of U for 
the FePt dimer (Figs. 7,8). It follows from these Figs, that with increasing U the position 
of the spin-down d x 2_ y 2 state changes significantly with respect to the Fermi energy (OeV), 
while that of other orbitals (or the position of the HOMO-LUMO gap) almost do not change, 
which means that their occupancies are not very much affected by U. The other important 
consequence of increasing U is an increase of the distance between the s- and d- states of 
the same spin, especially remarkable for the spin-up states. This separations indicates 
that one can neglect low-lying spin-up states for the s-orbitals for large U, except the d x 2_ y 2 
orbital (Figure 8). In other words, one can neglect the hopping processes between these and 
other states. In many cases, such a simplification may lead to a significant reduction of the 
computational time in the DMFT calculations. 

We have also anayzed the magnetic properties of four atom FePt systems. In particular, 
we have studied the U-dependence of the magnetization of the 3D and planar Fe 3 Pt by 
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using both the DFT+U and DFT+DMFT calculations (Figure 9). From Fig. 9 it follows that 
DFT+U leads to an increase of magnetization at large values of U for the 3D structures. This 
is not surprising, since in the case of 3D structure the correlation effects are more sensitive 
to the distance between atoms than in the 2D case in which the atoms are on average much 
further from each other even for small values of U. The reduction of magnetization due to 
dynamical effects, taken into account in the DMFT is even more dramatic than in the case 
of pure iron four-atom cluster (Fig. 4b). Indeed, when one iron atom is substituted by a 
platinum atom, the corresponding hopping parameter to this size increases significantly due 
to larger spatial extension of the corresponding orbitals, which decreases the average in-time 
occupancy of the orbitals, and hence decreases the magnetization. 

The chemical composition dependence of the magnetization for the four-atom clusters is 
presented in Fig. 10. In this case, the magnetization increases dramatically with increasing 
of the number of "more magnetic" Fe atoms as the DFT+U calculations show. The DMFT 
calculations show that that the magnetization is rather small and almost composition- 
independent when Pt atoms are present. This means that the hybridization effects are 
dominant even in the case of large Us. It is interesting to note that the DFT calculations 
for the Fe2-Pt2 cluster [13] give the magnetization 2/xb, which is much closer to our DFT+U 
result than to the DMFT one. This suggests that the energy levels calculated by the DFT 
approach are well-separated already at small U, and the hybridization is rather weak in this 
case. On the other hand, the dynamical effects may change dramatically the magnetization 
even in such a case of strongly-localized orbitals. 

V. POSSIBLE EXTENSIONS 

Since in most practical applications the clusters are supposed to be on a substrate, it 
would be very helpful to develop a DFT+DMFT formalism for such cases. One may think 
about this situation as a cluster coupled to a bath (substrate). Thus, for the lowest-order 
approximation one can solve the problem for the surface and than to take the substrate- 
cluster interaction into account by fixing the substrate energy levels. The inverse influence 
of the cluster on the substrate can be also taken into account as the next step, similarly to 
the approach proposed by Florens [16] . 

Another important case is the nonequilibrium behavior of nanosystems with strong 
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electron-electron correlations. In many technological applications, like electronic and mag- 
netic devices, usually fast switching of strong external fields takes place. Such a nonequilib- 
rium formulation of DMFT was already done in the case of extended systems [T2"]. In this 
case, one can study the time evolution of the local GF defined on the Kadanoff-Baym or 
Keldysh time contour. This problem is technically much more complicated than the time 
translationally-invariant case, since instead of solving one DMFT problem for each frequency 
in the equilibrium case, one needs to work with with complex GF and self-energy time ma- 
trices of a rather large size. The size of the GF matrix will be (number of atoms x number 
of orb itals) x (number of time points) instead of (number of atoms x number of orbitals) as 



in Eq. (14). Moreover, the number of time points must be twice the number of the real time 
points in the case of Keldysh contour, and the last number plus the number of (imaginary 
branch) points on the interval [0,1/T] for the temperature averaging in the case of Kadanoff- 
Baym contour. Typically, the number of the real time points in the nonequilibrium problem 
is around 1000 for the real branch and 10-100 for the imaginary one. Obviously, such large 
matrices can be barely inverted when one would like to solve the problem exactly, except 
the case of a few-atom cluster. In the case of nanosystems, one needs to use again the layer- 
by-layer solution proposed by Florens. Presently, we are working on developing numerical 
tools which will allow us to examine the magnetic characteristics of nanoalloys on substrates 
and also in presence of external fields. 



VI. CONCLUSION 



In this paper, we have provided and overview of the current status of the applications 
of combined DFT+DMFT approaches to study correlation effects in the case of molecules 
and nanostructures. In general, all the studies show that similar to the case of bulk and 2D 
systems, this approach is a promising method to take into account the correlation effects 
even in the case of small clusters and molecules. 

We have paid special attention to the case of nanomagnetism, for which we have proposed 
a DFT+DMFT approach and applied it to study the magnetic properties of small Fe and 
FePt clusters. It follows from our calculations that dynamical correlation effects lead to a 
significant decrease of the magnetization than DFT+U. There are two reasons for this: first, 
the time-dependence of the orbital occupancy taken into account in the DMFT approach 
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leads to a significant decrease of the magnetization relative to "the staggered" DFT+U case 
of "frozen" spins. Second, the frequency-dependence of the self-energy can lead to a shift 
of the energy levels with respect to the Fermi energy, and hence to a change of the level 
occupancy and magnetization. We have analyzed also the geometry, chemical composition 
and local Coulomb repulsion dependencies of the magnetization in the case of different 
clusters. In particular, we find that for a given system, to get the same magnetization from 
DFT+U and DMFT one needs to use a larger value of U in the latter, similar to results 
for the bulk. One of the reasons for this may be a stronger screening in the DMFT casejl]. 
The DFT+U and DFT+DMFT results are approximately on the same level of agreement 
with the experimental data. However, in the case of larger systems, one may expect that 
the DMFT results will much more accurate than the DFT+U. 

We have discussed a possible extension of the approach on two most important cases: 
particles and molecules on a substrate and nonequilibrium systems. The work in these two 
directions, including application to larger clusters, is in progress. 
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FIG. 4: DFT+U vs DMFT results for the magnetization of the Fe^ and Fe^ clusters at different 
values of U. 
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FIG. 5: DFT+U results for the magnetization of the Fe2 and FePt clusters at different values of 
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FIG. 6: DFT+U results for the positions of the s- and d-energy peaks in the DOS of the Fe2 and 
FePt clusters different values of U. The blue and red solid lines correspond spin up and spin down 
peaks of the Fe2 cluster. The dotted line correspond to the energy peaks for the FePt cluster. 
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FIG. 7: The density of states of different d-orbitals obtained within the DFT+U approach at 
U=0.1eV. 
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FIG. 8: The same as in Fig.7 at U=4.0eV. 
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FIG. 9: DFT+U vs DMFT results for the magnetization of the Fe^Pt clusters as a function of U 
in the case of 3D (a) and planar (b) geometries. 
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FIG. 10: DFT+U vs DMFT results for the magnetization of the four atom iron-platinum clusters 
as function of the number of iron atoms at U=4eV. 
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